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We compute the effect of an orbiting gas disc in promoting the coalescence of a central 
supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with 
explicit time dependence: we do not assume that the gas necessarily adopts a steady state or 
a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into 
a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the 
binary's gravitational radiation. We study the effects of star formation in the gas disc in a 
simple energy feedback framework. 

The disc spectrum differs in detail from that found before. In particular, tidal torques from 
the secondary black hole heat the edges of the gap, creating bright rims around the secondary. 
These rims do not in practice have uniform brightness either in azimuth or time, but can on 
average account for as much as 50% of the integrated light from the disc. This may lead 
to detectable high-photon-energy variability on the relatively long orbital timescale of the 
secondary black hole, and thus offer a prospective signature of a coalescing black hole binary. 

We also find that the disc can drive the binary to merger on a reasonable timescale only if 
its mass is at least comparable with that of the secondary black hole, and if the initial binary 
separation is relatively small, i.e. oq ^ 0.05 pc. Star formation complicates the merger further 
by removing mass from the disc. In the feedback model we consider, this sets an effective 
limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio 
is < 0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem 
for a significant class of mergers. 

Key words: accretion, accretion discs - black hole physics - galaxies: formation - cosmol- 
ogy: theory - instabilities - hydrodynamics 



1 INTRODUCTION 

In recent years, the process of shrinking a supermassive black 
hole (SMBH) binary by interaction with a circumbinary gaseous 
disc has been the topic of intense theoretical research. Initially 
this reflected attempts to overcome the 'last parsec' problem, 
i.e. the fact that dynamical friction with the stellar background 
is ineffect ive in shrinking the b i nary below s eparations smalle r 
than Ipc teegelman et alj 1 19801 : iMilosavlievic & Merritl 1200 lb . 
to the point where gravitational radiation can complete the co- 
alescence of the two holes. More recently there has been in- 
terest in finding an electromagnetic counterpart to the grav- 
itational wave emissio n expected during the final stages of 
black hole coalescence Jlvanov et all 19991: Armitage & NataraianI 



2002l:lMilosavlievic & Phinnevl2005l : lDotti et alj2006l : lLoebl2007F 
Cuadra et alj2009l) . 



In general most work on this subject has assumed driving by 
an accretion disc with constant mass inflow rate supplied from dis- 
tances far from the binary, and thus effectively assumed an infinite 



mass supply. However it is clear that in reality, where the gas is 
part of the galaxy merger producing the SMBH binary, the mass of 
gas ending up in the disc must be finite. Thus in general the disc 
does not settle to a steady state with a constant accretion rate, and 
its structure differs from a standard disc because it is affected by 
the tidal torque exerted by the binary. Moreover the rate or even 
the success of the shrinkage must depend on the total disc mass (as 
hinted by Cuadra et al, 2009). Here we investigate the problem by 
taking a disc of finite mass which is explicitly time-varying. We 
find significant differences from the results of assuming a steady- 
state accretion disc. 

We model the disc evolution in terms of a simple diffusion 
equation, including the effects of the tidal torques from the sec- 
ondary and calculating the binary orbit evolution self-consistently. 
Our initial condition is a finite disc mass concentrated at radii of 
the order of the initial binary separation. While still highly ide- 
alised, we regard this choice as more realistic than the assumption 
of a steady state, constant-M disc, into which a second black hole 
has somehow been inserted. 
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The paper is organised as follows. In Section 2 we describe 
the general features of binary shrinkage by a gas disc. In Section 
3 we describe our method for following both the binary and disc 
evolution. In Section 4 we describe the results. We constrain the 
disc mass required to bring the binary to coalescence and describe 
the expected appearance of such discs. In Section 5 we draw our 
conclusions. 



2 DISC ASSISTED BINARY SHRINKAGE 

The process of disc assisted binary shrinkage ('hardening') is for 
low mass ratios dynamically similar to the process of planetary mi- 
gration within a protostellar disc, and many studies have alluded to 
this analogy. For a low (but finite) mass companion, the tidal in- 
teraction between the binary and the gaseous disc produces a gap 
in the disc at the location of the secondary. The exchange of an- 
gular momentum across this gap is mediated through tides by the 
'satellite' (i.e. the secondary black hole) which behaves as a fluid 
element in the disc . In this case, the binary shrin ks on the disc's vis- 
cous timescale fy jArmitage & Nataraianll2002h . For a 10** Mq pri- 
mary black hole the viscous time at /? ~ 1 pc is a few times 10^ 
years, so the binary can shrink in a reasonable time before gravita- 
tional radiation emission takes over as the main cause of shrinkage. 

However things are different if the angular momentum associ- 
ated with the binary orbit is comparable to or larger than the local 
disc angular momentum (for a circular binary orbit, this is equiva- 
lent to requiring the secondary mass to be comparable to or larger 
than the local disc mass). Then the viscous torques in the disc are 
not able to redistribute efficiently the excess angular momentum 
transferred by the tides to the disc, and the binary shrinking slows 
down significantly. 

We define Mp and as the masses of the primary and sec- 
ondary blac k holes. In general the bi nary orbit evolution occurs on 
a timescale l lSver& Clarke, 1995. : .Ivanov et al.iri999i) : 



_ Mj(a) + M, 

^shrink , , , T 

Md(a) 



(1) 



which for Ms » M(j(a) can be significantly longer than the disc 
viscous time. In the equation above the local 'disc mass' is de- 
fined as Md(a) = 4nl,(a)a-, where E(a) is the surface density of 
the unperturbed disc at the binary separation a. In such a configu- 
ration, disc-assisted shrinking from parsec scales would take much 
longer than ~ 10' years and would effectively preclude binary co- 
alescence. Here the companion acts as a dam in the accretion flow. 
The disc inside the binary orbit is rapidly accreted, since it cannot 
be refilled from the outer circumbinary disc. Outside the binary or- 
bit, accreting gas builds up at the circumbinary disc's inner edge. 
This moves inwards on a timescale ~ fshrink, much longer than the 
viscous time. 

If we were now additionally to make the usual assumption that 
the disc is continuously fuelled from the outside, then the surface 
density at the inner edge of the circumbinary disc would increase 
until it wa s large enough to sp eed up the binary hardening, as de- 
scribed bv llvanov et al.l ( Il999h . The disc structure would be such 
that vL = M/2n = const, at large radii , while at sm aller radii 
vE oc typical of a 'decretion' disc jPringleiri99 ll) . where v 

is the disc viscosity. However, in reality, with a finite gas mass, the 
surface density at the inner edge of the c ircumbinary disc i s unable 
to grow steadily at the rate described bv llvanov et al.l ( ll99S^) . Evo- 
lution proceeds more slowly and a numerical solution is needed. 
We discuss these solutions in Section 3 below. The overall nor- 



malisation of the surface density differs from that predicted by 
Ilvanov et af] ( Il999l) . while the scaling with radius does follow the 
yS oc decretion-like solution. 

The qualitative appearance of the system as the binary evolves 
towards coalescence is clear. In a major merger two black holes 
of comparable ma ss are brought together at a distance < O.lpc 
dEscala et al]|2005l) . together with a finite amount of gas which we 
assume to settle rapidly to a fairly wide disc configuration. Part of 
the gas can accrete directly on to the primary black hole, producing 
a burst of AGN-like activity. Given the large amount of gas mass 
available, this burst is likely to be Eddington limited. As a result 
the disc around the primary black hole may be dispersed by the 
strong outflows arising from such a super-Eddington flow. Since 
the secondary hole acts as a dam for the gas outside the binary, 
the disc around the primary is not readily refilled. For protostel- 
lar accretion, the ext ent to whic h the dam is porous is still under 
discussion (e.g., Lubo w & D'An gelo 2006). Since AGN discs are 
typically much thinner than protostellar discs, the porosity of the 
dam for the AGN discs considered here is likely to be much re- 
duced, and we neglect this effect for the purposes of this paper. 
Thus after this initial episode, the system is characterised by a mas- 
sive circumbinary disc truncated at the binary separation, and pos- 
sibly by a low density circumprimary disc left over from the ini- 
tial Eddington-limited episode. The system appears bright at long 
wavelengths (optical/IR) qualit atively similar to the description of 
iMilosavlievic & PhinnevI i2005t) . 

At this point there are two possibilities: either (a) the circumbi- 
nary disc mass is too low, the coalescence fails (or takes inordi- 
nately long), and the two black holes are left orbiting at a fraction 
of a parsec distance, or (b) the disc mass is large enough to drive the 
binary to coalescence on a reasonable timescale. In this event the 
peak of the disc spectral energy distribution progressively moves to 
higher energies and the disc becomes more luminous, until gravi- 
tational radiation losses bring the black holes into coalescence. If 
the remnant circumprimary disc is still present, it will suddenly 
brighten as it is squeezed by the m erging binary, as qualitatively 
described bv lArmitage & Nataraiani 12002) . After coalescence the 
dam on the accretion flow is removed, so the inner disc hole is re- 
filled and the disc spectrum finally moves to high energies. 



3 TIME-DEPENDENT MODEL 

3.1 Coupled disc/SMBH binary evolution 

The evolution of an accretion disc in the presence of an embedded 
satellite can be described by a diffusive evolution model for the 
disc, including the tidal term arising from the secondary, i.e. 



dt RdR\ dR ) RdR^ ' 

where A = A/(nR)-, with A the specific tidal torque, and 



2 / \4 



2 \p 



R> a 



(2) 



(3) 



2 \P 

(We choose this formulation rather than the standard one with 
A, for later convenience.) In equation l[3j, CI ~ -sJGM^/W 
is the angular velocity at radius R, a is the radial position of 
the satellite, q = M^/Mp is the mass ratio between the sec- 
ondary and the primary black hole, and p = R - a. This sim- 



q' IRX* 
A= --^ - R<a. 
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plified form of the specific torque is commonly used in litera- 
ture (see, e.g., Armitage et al. 2002 : Armitage & Bonnelll I2OO2I : 
lArm itage & N ataraiaij|2002l : iLodato & Clarka l2004h . We smooth 
the torque term for R x a, where it would have a singularity 
(see equation 3 Tl. We use the same smoothing pres cription as in 
ISver & Clark37l995h and iLin & Papaloizoul j 19861) . i.e. for \R - 
a\ < max[//, /?h]i where H is the disc thickness and /?h is the size 
of the Hill sphere (Roche lobe) of the secondary. Finally we note 
that if the mass ratio q is large, the formalism above implies a non- 
vanishing torque at relatively large distances from the secondary. 
However, the tidal torque is actually prov ided by the cumulative 
effect of a series of Lindblad resonances jGoldreich & Tremaind 
and should therefore vanish beyond the outer and the inner 
Lindblad resonance, located at p = 2*-'^a. We therefore smoothly 
truncate the torque beyond these radii. 

The back reaction of the disc on the secondary orbital motion 
follows from angular momentum conservation, in the form 



df 



(4) 



-2nGM, 



f 



where the integral is taken over the whole disc surface and Q.^ is 
the angular velocity of the secondary. It follows that the behaviour 
of the system described here is determined by two dimensionless 
parameters. First, 

. QR- , 



-1 



(5) 



measures the relative strength of the second and first terms on the 
rhs of equation lO, while 



B 



Ani:(a)a^ 



(6) 



gives a measure of the magnitude of the rhs of equation l|4). 

The paramete r A has the simple physical interpretation 
jLin & Pap aloizou 1979h A = {l^jRf, where A is the gap width. 
To open a gap, the gravitational effect of the satellite has to over- 
come the pressure of the disc and the viscosity, which both oppose 
gap opening. We neglect the effects of pressure in this paper 

The second important parameter, S, represents the ratio of lo- 
cal disc mass to satellite mass. If the satellite is more massive than 
the disc (i.e. i? <k 1) its inertia makes migration slow. Conversely, 
if S » I, the secondary behaves like a fluid element of the disc and 
migrates on the local viscous timescale. 

3.2 Gravitational wave torques 

We also include in our calculation the torques arising from grav- 
itational wave radiation. For this purpose we add in Eq. l|4j the 
gravitational wave term causing the binary separation to evolve at 
a rate 



da 
df 



64 G'Ml 



5 c^a^ {\+qf 



(7) 



The corresponding timescale for a coalescence induced purely by 
gravitational wave emission is thus 

5 c^d" (l+qf 



64 G^Ml 



2 10' 



(8) 
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M„ 



O.lpc/ \lO**Mc 



a + q? 



-yr 



3.3 Energy dissipation from tides 

Most previous studies of satellite-disc interaction using a formal- 
ism like that above do not discuss in detail the extra energy dis- 
sipation introduced by the tidal term. As a consequence, the disc 
spectrum in the presence of a massive satellite is taken to be that 
of a standard accretion dis c truncated at a radius of the ord er of the 
binary separation (e.g., see lMilosavlievic & Phinnevll2005l) . 

In reality this point is subtle because the tidal interaction is 
not formally a dissipative term. However the formalism of Eqs. (2- 
4) enforces angular momentum conservation, and simultaneously 
assumes that the disc and secondary orbits remain circular. This 
implicitly assumes the presence of some dissipative term to damp 
out any eccentric motion. To compute this extra dissipation, we ig- 
nore the viscous term in equation ^ and follow the evolution of 
the system as driven purely by the tidal term. The total energy of 
the disc (kinetic plus potential) is: 



-27T 
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2R 



-mdR 



-TiGM, 
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and its time derivative is 
dt 



-nGM, 
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(9) 



(10) 



where in the last step we have integrated by parts. We now compute 
the change in the energy = -GMpMs/2a of the secondary. Using 
Equation (|4j we have 



d^ 
dt 



GMpM^da 
2cf dF 



-27rGMpfi, 



AI.dR. 



(11) 



We note that the energy changes of the disc and the satellite differ 
by a factor CI - Cls- Therefore the total energy loss due to the effect 
of tidal torques in Equations (2-4) is given by: 



dt 



: 27TGM, 



r 



(n, - Cl)AI.dR. 



(12) 



This energy loss is positive definite, as it should be, because both 
Q.S - ^ and change sign across the binary orbit. The local dissi- 
pation term is thus 



GM„ 



R 



-{n, - Q.)Ai.. 



(13) 



This term is only significant in a narrow range of radii at the disc 
edges, which are hotter than a simple truncated disc. We include 
this additional dissipation term in calculating the disc spectrum. 

This analysis assumes that the binary and disc orbits remain 
circular, and so would predict a bright circular disc rim. In reality 
the binary is likely to become somewhat eccentric. We discuss this 
point further in the concluding section of this paper 

3.4 Viscosity 

The final piece of physics needed in order to solve Eqs. (2- 
4) is a prescriptio n for the disc viscosity. We use the stan- 
dard a-prescription jShakura & Sunvaevll973h for a gas-pressure- 
dominated disc: 



(14) 
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where Cj = kT^/finiu is the sound speed at the disc midplane, is 
the midplane temperature, yu = 0.67 is the mean molecular weight 
and niH is the mass of the proton. The midplane temperature fol- 
lows from solving the vertical radiative flux equation for an opti- 
cally thick disc, assuming the opacity k to be either electron scat- 
tering or Kramers opacity, whichever the larger, i.e. 

Ivin' = a-sT*^ (15) 
o 

where (Tg is Stefan-Boltzmann constant, and T^g = Tc/t^^* is the 
effective temperature. The optical depth r is given by t = kL and 
the disc thickness H is calculated from the vertical hydrostatic bal- 
ance including both gas and radiation pressure. This setup, and in 
particular the requirement that the disc viscosity is proportional to 
the gas pressure, Eq. J14b . ensures t hat the dis c is thermally stable. 
We take a = 0.1 throughout cKing et al ]l2007l) . 



novae assu ming a Salpeter IMF from 0.1 Mq to IOOMq. This yields 
6 ~ 10"^ jThompson et al.l [2005h . However 6 could in princi- 
ple be lower or higher than this value. A higher value 6 ~ 0.01 
could be reached if the IMF of stars is top-heavy, as found for 
the star formation event in the central 0.5 pc of the Milky Way 
jNavakshin & Sunvaevll2005l : |Paumard et al.ll2006 l). A lower value 
of e might be appropriate if the stellar disc thickens quickly and 
adopts a geometrically thicker configuration than the gaseous disc 
(Nay akshin et al. 2007). Most of feedback energy would then be 
escaping from the disc rather than being deposited in the midplane. 

Another importan t effect is the increased pressure in the disc 
jThompson et al.ll2005h from winds and supernovae. However we 
limit our attention here to the energy effects of star formation only, 
and leave pressure effects to a future paper 



3.5 Gravitational instability and star formation 

The outer parts of AGN discs can be gravitationally unstable be- 
yond ~ 0.01-0.1 DC jKolv khalov & Sunvaev'l979':'Lodato & Bertiiil 
l2003l : lGoodmanfl2003l : iLodato 2007; King & Pringle 20o3). The 
gravitational stability parameter Q dToomrell 19641) is defined as 



Q 



(16) 



TCGI,' 

where kq is the epicyclic frequency. If 2 < 1 the disc is unstable. 
Further evolution of the disc is controlled by the ratio of the cooling 
timescale, fcooh to the dynamical timescale = fgooi t^. For < 1 the 
disc fragments into gravitationally bound objects jOamnniellioOll : 
[Rice et al. 2005). For longer cooling times, yS > 3, the instability 
saturates at a finite amplitude and the disc displays a long-lived spi- 
ral structure jLodato & Ric3l2004l2005l) . In protostellar accretion 
discs the parameter is generally much larger than unity ( iRafikovl 
I2OO5I) . However in AGN discs the opposite occurs (i.e./? « 1), and 
these discs are expected to fragment rather easily 1 Goodmaiill2003l : 
lNavakshidl2006l : iNavakshin etai]|2007l : lLevidl2007l) . 

The evolution of a fragmenting j8 <K 1 disc is uncertain. The 
prescription of a constant cooling time throughout the collapse of a 
fragmenting region is an oversimplification. If the collapse of self- 
gravitating fragments is not prevented, m ost of the mass quickly 
(almost dynamically) ends up in stars (e.g.. INavakshin et al.ll2007h 
and this might inhibit accretion altogether. Conversely, the energy 
input arising from the star formation process (either via protostel- 
lar accretion or from nuclear burning in the ne wly formed stars) 
could keep the disc close to marginal stability l lGoodmanll2003l : 
[Thompson et alj2005l : lNavakshinll2006h . 

In this paper we explore this latter possibility. In particular, at 
each timestep we compute the gravitational stability parameter Q 
(equationll6t for every radial bin in the disc. If 2 < 1, star forma- 
tion is allowed. We define the local star formation rate per unit area, 
Ssf . The appropriate value of this rate is found by requiring that the 
corresponding energy liberation rate is high enough to keep the disc 
in the state of marginal stability, assumed to be 2 = 1. The inter- 
nal heating of the disc is then dominated by the feedback from star 
formation, and hence the energy equation becomes 



(17) 



The free parameter of the model, e, is the efficiency of mass-to- 
energy conversion in star formation. The star formation rate per 
unit surface is added as a sink term on the rhs of Eq. l|2]l. 

A reasonable estimate of e is obtained by calculating the to- 
tal amount of energy liberated by nuclear reactions and super- 



4 TIME-DEPENDENT CALCULATIONS 

Here we describe the results of our simulations. We have run two 
sets, first without star formation (Section |4!2l l, and then taking feed- 
back from star formation into account (Section|4j3j- 



4.1 Initial conditions and disc setup 

We solve the set of equations (2)-(4) usi ng standard techniques 
(see, for example, lLodato&ClaA3l2004 . We use a logarithmi- 
cally spaced radial grid extending from the iimermost stable orbit 
at 6/?5, where = GM^jc^ is the primary gravitational radius, to 
an outer radius of 20ao, where ao is the initial binary separation. 
We typically use 100 mesh points. We use a zero torque boundary 
condition at both our innermost grid point (i.e. close to the primary 
black hole) and at our outermost one, at 20ao- The outer boundary 
condition does not really play a role since the disc never extends 
that far out. 

At each time during the simulation we compute the effective 
temperature Tsg(R) at the disc surface, including both the stan- 
dard viscous dissipation term and the additional tidally induced 
term described above. From T^g we calculate the disc spectral en- 
ergy distribution and the total luminosity by integrating a series of 
blackbody spectra with the appropriate radial dependent tempera- 
ture over the disc surface. 

All our calculations start with similar initial conditions. We 
take the primary black hole to have a mass of Mp = IO^Mq. We 
have run a number of simulations with different secondary mass, to- 
tal disc mass M^, and different initial binary separation Oq. The disc 
is initialised with a constant surface density over a narrow range 
of radii between 0.8ao and lag, where the density normalisation 
is chosen to give the required total disc mass. This is intended to 
give a simple realisation of the concept that both black hole and 
gas arrive at small radii as part of the same merger process. Al- 
though clearly an idealisation, we regard such initial configuration 
more plausible than the standard assumption that the black hole is 
somehow parachuted into an accretion disc that has already had 
time to relax into a steadily accreting state. We have considered the 
two initial separations ao = 0.01 pc and ao = 0.05 pc. During a 
gas-rich merger of two galaxies each hosting a black hole, the two 
holes can typically reach a distance of < 0. 1 pc by dynamical fric- 
tion against the s tellar and gaseous background l Escala et al .I2OO5I; 
iDotti et allbooel) . Thus our initial conditions are rather optimistic 
as to how far in a black hole can progress. Even so we find problems 
with merging on reasonable timescales even from such small initial 
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separations. This must imply stringent limits on tlie possibility of 
gas-driven coalescences. 



4.2 Discs without star formation 

In the first set of simulations we neglect the effects of star forma- 
tion and allow the disc to become unstable (2 < 1) without taking 
any further action. Although obviously unrealistic, these simula- 
tions are useful for comparing with those where star formation is 
included, and also to illustrate the state of the disc at decoupling, 
i.e. at the point where the binary coalescence starts to be driven by 
gravitational wave emission rather than by disc torques, as a func- 
tion of the main disc parameters. By insisting that no gas is removed 
by star formation we overestimate the degree to which the gas can 
facilitate a black hole merger. 

Our results shou l d be compared with those of 
iMilosavlievic & PhinnevI ( |2005[) , who consider the same gen- 
eral setup but do not follow the time-dependent evolution of the 
system. Instead they simply truncate the disc at the outer radius 
where it becomes gravitationally unstable. We therefore briefly 
summarise their main findings. They argue that, contrary to what 
happe ns when the mass ratio ratio ^ « I jArmitage & NataraianI 
l2002h . inner disc evolution for non-extreme mass ratios is much 
faster than the binary orbital evolution, so that after a short time 
the only gas present in the system is a circumbinary disc. We 
are interested in the state of the disc at decoupling, i.e. at the 
moment where binary evolution driven by gravitational wave 
losses becomes faster than the disc viscous evolution and the disc 
in ner edge cannot follow the seco ndary black hole as it spirals 
in. IMilosavlievic APhinne^ llool) assume the circumbinary disc 
at decoupling to be well described by a constant-M steady state 
solution, truncated at the inner edge. Using the requirement that 
the disc is thin at the irmer edge, they argue that M < MEdd> where 
Mem = i'Edd/'/c^ is the Eddington accretion rate evaluated at the 
innermost stable orbit, Lndd is the Eddington luminosity and 77 is 
the accretion efficiency. The resulting spectral energy distribution 
is therefore approximately given by a standard multi-colour 
blackbody spectrum, but lacking the high-energy component and 
therefore peaking at optical-IR wavelengths. After decoupling 
the binary torques are removed and the disc is free to flow to the 
bottom of the potential well, producing in a bright flare, and thus 
recovering the high-energy part of the standard disc spectrum. As 
we shall see, this is not the complete picture. 

4.2.1 Initial binary separation equal to 0.01 pc 

We consider first the case where the initial binary separation is aq = 
0.01 pc, although the general behaviour of the system is similar 
for all the simulations. At this distance from a 10* Mg primary, the 
timescale for a merger induced purely by gravitational radiation is 
fg„ w 2 X 10'' q-' yrs. 

The luminosity as a function of time for all our simulations is 
shown in Fig.[T] where the left panel refers to q = 0.3, the middle 
panel refers to g = 0.1 and the right panel refers to ^ = O.OI. The 
various lines indicate different values for the total disc mass 
(see caption for details). We consider a range of cases varying from 
a disc mass comparable to the secondary mass M^/M, = 1 to one 
where it is much smaller M^/Mj = 0.01. When the disc mass is 
much larger than the secondary mass, the secondary black hole is 
simply swept in to the central black hole by the accretion flow. 

We start with roughly half of the gas inside the orbit of the sec- 
ondary black hole and half outside. Initially, the matter inside the 
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Figure 3. Surface density of the outer disc for the case y = I, /Ms = 1 . 
The lines refer to / = 0.5 Myr (solid line), 0.75 Myr (short-dashed line), 
IMyr (long-dashed line), 1.5 Myr (dot-dashed line) and 1.6 Myr (dotted 
line). The last two cases refer to epochs when the binary orbit accelerates 
because of gravitational wave emission (cf. Fig. |2), and the torque on the 
inner disc is removed. 



binary orbit accretes on to the primary black hole, resulting in a rel- 
atively bright AGN, where the luminosity obviously scales with the 
total disc mass. For the choice of parameters made here the peak lu- 
minosity is at most of order the Eddington luminosity for our most 
massive disc. The inner disc evolves faster then the binary orbit and 
the surface density and the disc luminosity drops rapidly. However, 
a low density disc is still present at decoupling in all of our simula- 
tions. Clearly, the amount of mass left over in the inner disc at de- 
coupling is strongly model-dependent. However, the presence of a 
small, but non-negligible, mass in the inner disc might have impor- 
tant observational consequences. As explained above, even if such 
a low density disc does not provide a significant contribution to the 
SED, it can be squeezed when angular momentum loss by gravita- 
tional radiation becomes dominant, producing a sudden burst of lu- 
minosity. After this, the outer disc can finally flow to the innermost 
region and provide a longer-lasting AGN-like appearance, even 
though the post-merger luminosity is generally sub-Eddington. 

Fig.|2]shows the evolution of the binary separation with time. 
The initial evolution is rather slow, being driven mostly by the disc 
torques. It suddenly accelerates when the gravitational radiation in- 
duced evolution takes over and the outer disc decouples from the 
binary evolution. We see from Fig. |2] that the radius at which the 
disc decouples increases with decreasing disc mass. A similar re- 
sult was also found by Milosavlievic & Phinnev (2005), who pa- 
rameterised the disc mass terms of the mass flux at the disc inner 
edge (see below). Note that for disc masses smaller than roughly 
O.I of the secondary mass, the merger time becomes comparable 
with the gravitational radiation timescale, indicating that the disc 
gives only a limited contribution towards accelerating the merger. 

Figure [3] shows a series of snapshots of the surface density in 
the outer disc for the case where ^ = 0. 1 and Md/M, = 1 (which 
corresponds to the solid line in the middle panels of Figs.[T]and|2]l. 
The lines refer to t = 0.5 Myr (solid line), 0.75 Myr (short-dashed 
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Figure 1. Luminosity evolution of circumbinary discs for an initial separation of oq = 0.01 pc. The secondary/primary mass ratio is Ms/Mp = q = 0.3 (left 
panels), 0.1 (middle panel) and 0.01 (right panel). For each panel, the various lines refer to different disc masses: Mj/Mj = 1 (solid lines), 0.5 (dotted line), 
0.1 (short-dashed line) and 0.01 (long-dashed line). 




Figure 2. Evolution of the binary separation a for the initial value of c/q = 0.01 pc. The three plots have the same parameters as the coiresponding plots in Fig. 
[T] Similarly, the line styles are the same as in Fig.[T] 



line), IMyr (long-dashed line), 1.5 Myr (dot-dashed line) and 1.6 
Myr (dotted line). This clearly shows that in the situation consid- 
ered in this paper the surf ace density at the inner edge does not 
grow with time (unlike the llvanov et al ] | 19991 solution). The even- 
tual growth of E at the inner edge only occurs at late times, when 
the binary orbit accelerates due to gravitational wave emission (cf. 
Fig.|2j and the tidal torque is effectively removed. 

Figure |4] shows a number of properties of the disc at decou- 
pling, which is the moment at which the binary separation evolution 
becomes faster than the disc viscous time close to the binary orbit. 
The top panels show the disc surface density profile (left) and the 
profile of the local mass flux in the disc M = InvL in units of the 
Eddington rate at the innermost stable orbit M^^, where we assume 
mass to energy conversion efficiency of r] = Q.\ (we use th is scaling 
for ea se of comparison with the results of lMilos avlievic & PhinnevI 
I2OO5I) . The bottom plots show the aspect ratio HjR (left) and the ef- 
fective temperature profile (right). These plots refer Xo q = OA and 
the three lines correspond to to M^/M^ = 1 (solid line), 0.5 (dotted 
line) and 0.1 (dashed line). 

There are several interesting features in these plots. We first 
discuss the radial profiles of M. The main thing to note is that the 
outer circumbinary disc does not have a constant M, as a stan- 
dard steady accretion disc model would require. The tidal torque 
at the inner edge of the circumbinary disc changes its overall struc- 



ture, such that if a steady state is reached the mass flux scales as 
M oc This kind of profile is typical of visc ous disc solu- 

tions in the presence of a torque at the inner edge ( |Pringlelll99ll 
ISver & Clark3ll995l : llvanov et alj|l999t) . and as we will discuss 
below it affects the shape of the spectral energy distribution of 
the disc. We stress that the choice of scaling M to the inner Ed- 
dington rate is motivated by convenience, and the mass flux can 
significantly exceed the reference value (the Eddington value for 
the innermost stable orbit) without the luminosity becoming super- 
Eddington, since the disc is truncated at a radius much lar ger than 
the innermost stable orbit. iMilosavlievic" & PhinnevI l l200^ still ar- 
gue that the limit M/M^id < 1 is required in order to keep the disc 
thin at its inner edge, at least for relatively large q. Although our 
most massive disc does have Ml Mem ~ 1 at its inner edge, we 
find that the disc aspect ratio HjR is <k 1 (see bottom left panel), 
therefore in principle allowing a much larger inner mass flux with- 
out violating the thin disc requirement. The straight (red) line in the 
top right panel of Fig.|4]shows the scaling between the inner edge of 
the circumbinary disc at decoupling and the mass flux predicted by 
Milosavlievic & Phinnev (2005), where the normalisation has been 
chosen to match our results. Although our simulations do repro- 
duce this scaling quite well, our results imply a binary separation 
at decoupling larger by a factor two than lMilosavlievic & PhirmevI 
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Figure 4. Properties of the disc at decoupling, i.e. wlien the binary evolution via gravitational wave radiation becomes faster than viscous evolution at the 
secondary's orbit. The initial separation was ciq = 0.01 pc. The secondary/primary mass ratio is ^ = 0.1. The different hnes refer to different values of the 
initial disc mass, with the same notation as Fig.[T] Shown here is: the disc surface density profile in CGS units (upper left), the mass flux through the disc in 
units of the Eddington rate at the innermost stable orbit (upper right), the aspect ratio H/R (lower left) and the effective temperature in Kelvin (lower right). 



( l2005h . At decoupling we have a = 335, 520, 1400i?g, for the three 
cases shown in Fig.lH i.e. Mi/Ms =1, 0.5 and 0.1, respectively. 

The top left panel in Fig. |4] shows the surface density profile 
at decoupling. Here we note that although the inner circumprimary 
disc has been significantly depleted in mass and its surface density 
is thus low, it has still not disappeared. As emphasised above, even a 
small amount of mass in the circumprimary disc produces a bright 
burst as it is squeezed by the final evolution of the binary orbit, 
producing the sudden peak in the light curves plotted in Fig.[T] The 
bottom right panel in Fig.|4]shows the effective temperature profile 
at decoupling. Here the most significant feature is the presence of 
two bright spots at the edges of the gap around the secondary. This 
results from the extra dissipation term discussed in Section 3.3, and 
provides a significant, and previously neglected, contribution to the 
overall spectral energy distribution. 

The resulting spectral energy distribution (SED) for the three 
cases of Fig.|4]are shown in the top panel of Fig.|5] These are calcu- 
lated as multicolour blackbody spectra for the appropriate effective 
ternperature pr ofiles. The general shape is similar to that predicted 
bvlMilosavliev ic & Phinn ev (2005), although there are significant 
differences. The disc is very bright at optical/IR wavelengths and 
lacks the high energy emission from the inner disc. The inner disc 
emission is generally negligible because of the low density of this 
disc, although some high energy emission from the inner disc can 
be seen for the lowest disc mass (dashed line). The bottom plot in 
Fig. [5] shows the contribution to the total SED of the various disc 
components, for M^/M^ = 1. The short-dashed line shows the con- 
tribution of the outer circumbinary disc only, without including the 
extra dissipation term at the gap edges. The long-dashed line shows 
the spectrum arising from the whole disc, but neglecting again the 
dissipation at the gap edges. This emphasises the contribution of the 
inner disc. The solid line shows the total SED, including the emis- 
sion from the gap edges, which in turn result in a significant contri- 
bution at relatively higher energies. The red line, for comparison. 



shows the expected SED of a truncated d isc with constant M. This 
is the shape of the SED as predicted by Milosavlievic & PhinnevI 
( l2005h . As one can see, the long-wavelength part of the SED of a 
constant M disc is shallower than the one we calculate. The steep- 
ening of the SED in t he presence of a massive satellite had been 
discussed by ISver & Clarke ( 1995). Then we have vFy oc y ' , 
rather than the value of vFy oc v"*'^ found for a steady accretion 
disc. 

Note that the frequency range of the hot edge emission is com- 
parable to that of the inner disc. Indeed Fig. |4] shows that the ef- 
fective temperature of the edge is similar to that of the inner disc. 
However, the area covered by the hot edges is much larger than the 
inner disc and provides a much larger luminosity. The total lumi- 
nosity of the hot edges is roughly equal to half of the bolometric 
luminosity. 

After the black holes coalesce and the tidal torque at the in- 
ner edge is removed, the gas will eventually freely flow to the in- 
ner disc, providing a luminous high energy component to the SED, 
which would gradually flatten to reach a standard constant-M so- 
lution in the inner disc. 



4.2.2 Initial binary separation equal to 0.05 pc 

For a larger initial binary separation of aq = 0.05 pc, the grav- 
itational radiation timescale is far longer and would not produce 
coalescence in a Hubble time without any further binary hardening 
mechanism. Here again we take the primary black hole mass to be 
10** Mq. The mass ratio between the secondary and the primary is 
^ = 0. 1 and we consider, as in the previous section, different values 
for the disc/secondary mass ratio. Fig.|6]shows the evolution of the 
binary separation for Mil Ms = 1 (solid line), 0.5 (dotted line) and 
0.1 (dashed line). Here again, we see the same trend as observed 
in for the smaller initial separation Aq = 0.01 pc. We note that disc 
torques are able to induce coalescence within a Hubble time only 




Figure 5. Top: SED of the disc (in CGS units) at decoupling for various disc 
masses: Mi/M^ =1 (solid line), 0.5 (dotted line) and 0.1 (dashed line). The 
initial separation was ao = 0.01 pc. Bottom: Contributions to the overall 
SED for disc mass M^/Ms = 1 shown in the top plot. The solid line indicates 
the full spectrum; the short-dashed line shows the contiibution of only the 
outer disc, without including the extra tidal dissipation term; and the long- 
dashed line includes the emission of the inner disc but still does not include 
the hot gap edges. Comparing the long-dashed and solid lines thus shows 
the importance of the hot-gap emission at high radiation frequencies. For 
reference, the red line shows the SED of a standard constant-M accretion 
disc: its spectrum is shallower than the one found here, which more closely 
approximates the structure of a decretion disc. 



when the disc mass is comparable to the secondary black hole. At 
the end of our simulations (~ 1 Gyr) the binary had merged only 
for Mi = Mj. The final state of the disc at decoupling follows the 
same trend as before. 



4.3 The eflfects of star formation 

From the analysis of the results of the previous section, we may 
draw the following general conclusion. Even starting from an initial 
separation of ciq = 0.05 pc, a circumbinary disc is able to bring the 
SMBH binary to coalescence in a reasonable time only if the disc 
mass is of the order of, or slightly smaller than, the secondary black 
hole mass. However, in the previous section we did not include 
the effects of gravitational instability and simply allowed the disc 
to have a value of Q smaller than one in its outer parts, beyond 
~ 0.05 pc. Here, we include the effects of gravitational instability 
by requiring that, whenever the disc is unstable, it is kept in a state 
of marginal stability by the feedback induced by star formation. 
We thus remove from the calculation the amount of mass needed 
to be converted into stars to provide such feedback, as detailed in 
Section [331 The effect of removing gas from the disc is to reduce 
the efficacy of the disc in merging the black holes. 

We first consider the case where the secondary/primary black 
hole mass ratio is q = 0.1, the initial mass ratio between the disc 
and the secondary is M^/Mg = 0.5 and the initial binary separa- 
tion is 0.01 pc. We also set the star formation feedback efficiency 
to e = 10"^. In Fig. [7] we show the resulting evolution of the bi- 
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Figure 6. Evolution of the binary separation with time when the initial 
separation is taken to have the larger value of oq = 0.05 pc. The primary 
mass here is lO^M©, the secondary/primary mass ratio is ^ = 0.1 and the 
disc/secondary mass ratio is Mi/Ms = 1 (solid line), 0.5 (dotted line) and 
0.1 (dashed line). Only the more massive discs produce a merger within 1 
Gyr. 



nary separation (solid line) and that in the corresponding calcu- 
lation without allowing for star formation (dashed line). Star for- 
mation significantly reduces the ability of the disc to induce the 
merger. This is because star formation depletes the gas reservoir in 
the outer parts of the disc, reducing its surface density, and so sig- 
nificantly reduces the disc density just outside the secondary orbit. 

Fig. [8] shows the evolution of the total star formation rate 
within the disc. We can see that over the first 10^ years the disc 
undergoes significant star formation. The total mass converted into 
stars is of the order of » 2 10* Mq, that is roughly half of the ini- 
tial disc mass. Finally, we show in Fig.|9]the main properties of the 
disc after a 6 Myr from the beginning of the simulation, after most 
of the star formation has died out but before the secondary has mi- 
grated significantly. The top panel shows the disc surface density, 
the middle panel shows the aspect ratio H/R and the bottom panel 
shows the value of the stability parameter Q. We can thus see that 
our star formation prescription indeed causes the disc to hover near 
marginal stability in its outer parts. 

Summarising our findings for the simple energy feedback star 
formation model, we note the following. In this model, the evolu- 
tion of the system does not depend significantly on the initial disc 
mass (except when it is so low as to be gravitationally stable at all 
radii, and would not be able to drive the merger). Thus different 
initial disc masses only imply differing amounts of star formation 
rate at early times, and the surface density after a few million years 
is very similar in all cases. 

The effect of changing the efficiency of star formation feed- 
back 6 is not trivial. Naively one would expect to simply have more 
star formation as one decreases its value (cf. eq. JI7I(). thus lead- 
ing to a smaller disc density E. However, the reduced E requires a 
smaller temperature in order to be gravitationally stable, which in 
turn implies a much smaller cooling rate and a lower star formation 
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Figure 7. Evolution of the binary separation with time in calculations with 
(solid curve, see ^4.3) and without (dashed, see ^4.2t star formation in the 
disc. The latter case is the one shown in Figures|2]and|4] The primary mass 
is 10^ Mq, the secondary/primary mass ratio is ^ = 0.1, the disc/secondary 
mass ratio is Mj/Mj = 0.5 and the initial separation is uq = 0.01 pc. Note 
that removal of gas by star formation significantly slows down binary merg- 
ing. 



rate. We have run different simulations with e = W^^ and e = lO""* 
and found only found small differences. 

Finally, in Fig. [TO] the solid line shows the evolution of the 
binary orbit in the presence of star formation when the initial sep- 
aration is 0.05 pc, the secondary/primary mass ratio is ^ = 0. 1 and 
the initial disc mass is equal to the secondary mass. This should 
be compared with the equivalent run where star formation is not 
included (dotted line). Clearly, while in principle there is initially 
enough gas mass in the disc to bring the holes to coalescence, the 
effect of star formation depletes the disc severely and the amount of 
gas left over is not enough to lead to coalescence within a Hubble 
time. 



5 DISCUSSION AND CONCLUSIONS 

We consider the merging of two black holes in the centre of a 
galaxy driven by a gaseous disc. We make the assumption that 
the secondary black hole (mass Ms) arrives close to the already 
present p rimary (mass M„) as a result of a galactic merger of some 
kind (cf. lKing& PringI bringing with it a certain amount 

of gas (mass Mj). This contrasts with the usual assumption that the 
merger process manages to parachute the secondary black hole into 
an already existing steady state accretion disc of infinite extent and 
mass. 

If the amount of gas is large (Mj » M,) and does not mainly 
turn into stars then the secondary black hole is simply swept into 
the primary by the resulting accretion flow. Thus we consider what 
happens when the amount of gas is comparable to, or smaller than, 
the secondary mass. We first ignore the possibility of star forma- 
tion within the accretion disc flow (Section|4j2]( in order to demon- 
strate the time-dependent properties of the interactions between 
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Figure 8. Total star formation rate for the simulation shown also in Fig. [7] 
Over the entire course of the simulation the total amount of gas turned into 
stars is » 2 IO'^Mq. 
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Figure 9. Disc properties after 6 Myr At this time star formation has con- 
siderably depleted the disc mass and little disc evolution has taken place. 
Because of this the results are relatively independent of initial disc mass. 
The initial separation was oq = 0.01 pc, the secondary/primary black hole 
mass ratio is q = 0.1, the disc mass is Mj/Mj = 0.5 and the star formation 
feedback efficiency is e = 10"^. Top panel: disc suiface density in cgs units. 
Middle panel: aspect ratio H/R. We see that typically H/R a 10"'. Bottom 
panel: gravitational stability parameter Q. 
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Figure 10. Same as Figure |7] but for the initial separation uq = 0.05 pc 
and tlie disc/secondary mass ratio is M^/Ms = 1. As before, solid curve 
shows the case with star formation included, whereas the dotted shows that 
where star formation ignored. Star formation severely depletes the disc, and 
the remaining disc mass is not large enough to induce a black hole merger 
within a Hubble time. 



disc and secondary black hole. We then, more realistically (Sec- 
tion l4.3l ). take account of formation of stars from the disc gas when 
it becomes self-gravitating, and take account of the feedback from 
those stars on the gas in the disc. 

We have computed the properties and the SED of the disc 
at and around the time of decoupling (i.e. when gravitational 
radiation takes over as the dominant merging effect) in cases 
where the merger occurs. Our findings differ from earlier results 
jMilosavlievic & PhinnevlllOOSi) for three main reasons. 

(a) The disc is not in a steady state with uniform M, but closer 
to a decretion disc with M oc i?"'/^. This makes the SED steeper 
(vFy oc y", rather than the steady accretion disc value of vFy oc 
y'-3). 

(b) The region inside the secondary orbit might not be empty, 
and the small gas density can still lead to a strong burst of luminos- 
ity durin g the final stages of the merge r, as predicted for low mass 
ratios bv lArmitage & NataraianI j2002h . 

(c) The enhanced dissipation resulting from the interaction be- 
tween the disc and the black hole heats the disc edges either side 
of the secondary significantly. This is a new and general result, im- 
plying an extra high photon energy contribution to the SED. 

The existence of this harder e mission weakens the prediction 
of iMilosavlievic & PhinnevI bOOSh that the presence of a binary 
can be easily inferred by a shift of the peak emission to low pho- 
ton energies. More positively, it offers a new way of inferring the 
presence of a binary through a com bination of this resu lt with pos- 
sible detectable periodic variability l lHaiman et al.l2009l) . Although 
in this paper we have assumed that the secondary and disc orbits are 
exactly circular, in reality the system may be driven to somewhat 
eccentric orbits. For disc masses large enough to shrink the binary 
significantly (as assumed here), dissipation in the disc is likely to 
keep the eccentricity low, and so the general picture remains similar 
to that we have described. However the dissipation at the inner edge 



of the disc will not be uniform in azimuth or in time. This implies 
variability at photon energies which one would normally associate 
with the innermost disc region around a single black hole. How- 
ever the typical timescale of this disc-gap emission is not that of 
the innermost orbits around the primary black hole, but instead the 
much longer orbital timescale of the secondary. Thus, the detection 
at high photon energies of variability on the long timescales typical 
of the outer disc would be a clear signature of the presence of a 
black hole binary. 

Evidently detailed modelling of such effects requires the use 
of 2D (or even 3D) disc models, as it depends on the effectiveness 
of disc waves generated by tides to carry away the excess tidal en- 
ergy. 

Circumbinary discs are gravitationally unstable at a O.lpc, 
just as their cousins around single SMBH are. It is not yet clear 
if star formation "catastrophe" can be avoided due to the action of 
thermal or momentum feedback from m assive stars in the disc (e.g., 
lGoodmanll2003l : [Thompson etd]|2005 ^ . In the case of the thermal 
feedback, as shown here, gas is converted into stars too quickly 
compared to the binary evolution timescales. This has a severe im- 
pact on the ability of the disc to induce a black hole merger in a rea- 
sonable time. We find that the rapid effect of star formation sets an 
upper limit to the effective mass of a gaseous disc. Because of this, 
mergers within a Gigayear require initial separations closer than 
ao = 0.01 pc rather than 0.05 pc, implying that gas driven mergers 
of this sort are unable to solve the 'final parsec' problem. To avoid 
the disc becoming self-gravitating the disc must have Q > I and 
so must have a mass at most M^/Mj, « H/R. We also require that 
the disc mass be at least comparable to the secondary mass to pro- 
duce a merger. Our modelling of star formation leads to disc aspect 
ratios H/R of order of a few times 10"' and so we conclude that 
disc assisted mergers only work for mass ratios q < 0.001 if only 
thermal feedback is included. Whether such low mass secondaries 
would be able to produce a ti ght pair in the first place is a matter of 
debate jCallegari et"Zll2009l) . 

Our treatment of star formation is very simplistic by necessity. 
The inclusion of other effects neglected here, such as momentum 
feedback and dynamical support from the newly born stars embed- 
ded in the disc, might stabilse the disc at higher surface densities. 
Additionally, the process of fragmentation might result in the for- 
mation of clouds and clumps in the disc, rather than stars, whose 
dynamics is still unclear. Both effects would results in a relatively 
more massive gaseous disc at ~ 0.1 pc. Finally, dynamical fric- 
tion from the newly born stars (not included here) might help in 
bringing the holes together. While it is hard to make a precise esti- 
mate, we believe our calculations with and without star formation 
(see Figures[7]and|10l) bracket the range of possible timescales for 
disc-induced binary mergers. 

Nevertheless it seems quite likely that the 'final parsec' prob- 
lem is difficult to overcome in all cases, particularly for non- 
extreme mass ratios. This would imply the relative ubiquity of sec- 
ondary black holes in galactic nuclei. 
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